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Stochastic processes associated with traveling wave sohitions of the sine-Gordon equation 
are presented. The structure of the forward Kolmogorov equation as a conservation law is 
essential in the construction of the stochastic process as well as the traveling wave structure. 
The derived stochastic processes are analyzed numerically. An interpretation of the behaviors 
of the stochastic processes is given in terms of the equation of motion. 
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CN ■ 1. Introduction 
> , 

Tlj" ' Research on stochastic processes has received much attention in mathematical physics 

T^lj- . recently. For example, the stochastic Lowner equation attracts much interest regarding its 

connection to the conformal field theory,^ as well as its original meaning as a stochastic version 
ON ■ 

' of the Lowner equation.'^ The stochastic cellular automaton and asymmetric simple exclusion 
00 ! 

I process (ASEP), among others, are also attracting considerable interest. The ASEP describes 



a non-equilibrium process of interacting particles, whose steady state has been exactly ob- 
^ ■ tained."^ Even the current at equilibrium and the corresponding phase diagram were exactly 

■ presented for the ASEP.^ Besides the ASEP, several stochastic cellular automata are invented 

for the successful analysis of traffic flows. ^ 

The stochastic process associated with the Burgers equation is called Burgers' process,^ 
which is applied to studies of turbulence.* This process is a typical example of nonlinear 
stochastic processes, in other words, stochastic processes associated with nonlinear equations. 
The time evolution of the density function of Burgers' process, for instance, is governed by 
the Burgers equation Uu — 2uux + Uxx = 0, which is indeed nonlinear. Such nonlinear stochas- 
tic processes are expected to be an effective approach to analyzing phenomena, where both 
nonlinear and stochastic effects are simultaneously predominant. 

Among many kinds of nonlinear equations that describe various physical systems, we have 
a series of integrable equations. Expecting future applications to the analysis of phenomena 
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related to integrable systems, we shall present new nonlinear stochastic processes associated 
with the sine-Gordon (SG) equation 

I?!"!* = 'Ti^ sin ^, m: constant, (1-1) 

along the same line as the Burgers process. As is well known, the SG equation is an integrable 
nonlinear equation that is suitable for describing the dynamics of connected pendulums under 
gravity,^ the motions of the dislocation of one-dimensional materials, the propagation of 
optical waves in one-dimensional space,^^ and so forth. The SG equation possesses a kink 
solution that has localized structure, aside from periodic solutions. In this paper, we shall 
consider periodic and one-kink solutions having traveling wave structures, and we shall present 
stochastic processes associated with these solutions. We shall also perform numerical analyses 
of the derived evolution equations under suitable initial conditions in order to study their 
behaviors. The results and their explanation given in the later sections are expected to describe 
typical features of stochastic motion associated with integrable equations. 

This paper is organized as follows. In the next section, the stochastic differential equation 
is presented. A detailed description of the derived stochastic processes is given in §3. In §4, 
numerical calculations of the stochastic equations for some kinds of traveling waves are pre- 
sented, and interpretations of the results of the numerical analyses are shown in §5. The final 
section is devoted to conclusions, including discussions on possible experimental realization. 

2. Forv^ard Kolmogorov Equation as a Conservation Law of the Sine-Gordon 
Equation 

We shall consider a one-dimensional Ito diffusion, which is a scalar stochastic variable 
Xf obeying the stochastic equation 

dXt = b{Xt, t) dt + a{Xt, t) dBt, (2.1) 

where the real-valued functions b and a are called the drift and diffusion coefficients, respec- 
tively. The one-dimensional Brownian motion is denoted Bf. The density function t) obeys 
the forward Kolmogorov equation 

^ = A*p, (2.2a) 

where A* is defined by 

A* = -^ibix,t)-] + l^[a\x,t).]. (2.2b) 

The operator A* is called the conjugate operator of the generator A of eq. (2.1). Equation 
(2.2a) can be interpreted as a conservation law ofp{x,t), with the flux —bp+ ((T^p)^/2. 
Among the conservation laws of the SG equation, we shall consider 

{ct)l), = i-2m^cosct>),. (2.3) 
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Equations (2.2a) and (2.3) have the same form when we define p = N(p'^, where A'^ is the 
normahzation factor. Thus, we can read 

-Hi + = -2m2 cos ct> + A{t), (2.4) 

where A{t) is an arbitrary function. When is a travehng wave solution, 

(f) = (p{x — vt), v: constant, (2-5) 

the drift and diffusion coefficients are separated into simple forms. Since the derivatives of the 
variable are given as 

(l>x = 4>', <t>t = -v(t)', (2.6) 

and (j)xt = we find that —v<p" = m^sin^ from eq. (1.1). Multiplying both sides by (j)' 

and integrating by x, we have 

(j)' =4>l = -^(cos - C), C: constant. (2.7) 
Thus, we can eliminate coscf) from (2.4) and derive 

"-^^ ^{a^4>l)x + A{t) = 2m'^C + v(l)l. (2.8) 



2 



Introducing a new function F = a (f)^, we can rewrite eq. (2.8) as 

Because b and a should be bounded for all x values even at the zeros of (px, we assume that 
F and A are given by 

F{x,t) = A = Am'^C. 

Then the coefficients of eq. (2.2a) are determined by 

6=^+f<Axx</'r^ ^' = ;^'^r'' "^2. (2.9) 

The above restriction to n prevents b and a from diverging. Hence, the stochastic differential 
equation eq. (2.1) becomes 



dXt =(v + |0..0r') dt + dBt, (2.10) 

where v is the velocity of the traveling wave, n > 2, and a are constants. Equation (2.10) 
can be considered to describe the stochastic motion of a particle moving under the influence 
of a traveling wave solution of the SG equation. When the function cp^ approaches zero ex- 
ponentially, such as the kink solution, eq. (2.10) is also well-defined for 1 < < 2, since 
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3. Traveling Wave Solutions of the Sine-Gordon Equation and the Coefficient of 
Associated Stochastic Differential Equation 

Wc shall consider the coefficients in eq. (2.10) for several types of solutions (f) of the SG 
equation (1.1). In eq. (2.7), the velocity v is allowed only negative values for C > 1 and 
only positive values for C < — 1, whereas both positive and negative v values are allowed for 
\C\ < 1. Among the traveling wave sohitions, we choose those that have negative velocities: 



^(z) = 2 arctan 



2 arctan < 



4 arctan 



C-1 
C + 1 
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(3.1a) 



(3.1b) 



(3.1c) 



The functions cq{z, k) and tn(2;, k) are Jacobian elliptic functions with a modulus k. We note 
that the solution is periodic for C > 1 and |C| < 1, whereas the solution forms a solitary wave 
for C = 1. In the case associated with periodic solutions, their supports need to be compact 
for these cases or a periodic boundary condition for x is required, since the normalization 
factor of the density function should be finite. The solution eq. (3.1c) is the kink solution, 
and we can consider the stochastic equation on the whole space. Hereafter, the parameter n 
is chosen to be either n = 1 or 5. In the former case, its stochastic effect is spatially uniform. 
In the latter case, on the other hand, its stochastic effect depends linearly on the density 
function of Xf. 

In Figs. 1 and 2, we show two examples of the drift coefficients a in the case of n = 1, 
and we display no diffusion coefficients because they are constants. Figure 1 shows a typical 
configuration of 6 at f = for the periodic solution eq. (3.1a) with C > 1. The parameter 
C cannot have a value at |C| < 1, because periodical zeros of (jjx in the denominator of b 
yield singularities in b, and the stochastic process is not well-defined. Since b is defined by eq. 
(2.9), it is periodic owing to the periodicity of 4>. In the case associated with the kink solution, 
C = 1, the figure of b is given in Fig. 2. The coefficient b has a localized structure around the 
center of the kink. 

In Figs. 3 and 4, we shall show two typical examples of the drift and diffusion coefficients 
in the case of n = 5. Figure 3 shows the coefficients for periodic solutions C > 1. Both 6 and a 
have the same period, which is half of the period of the solution (f). On the other hand, in the 
case of |C| < 1, the appearances of b and a are similar to those shown in the case of C > 1, 
but their period is the same as that of <p. 

In Fig. 4, we show the coefficients of eq. (2.10) for the kink solution (3.1c). Both of the 
coefficients have localized structures around the center of the kink. Stochastic behavior of the 
variable Xf appears only in the vicinity of the kink, and the drift coefficient changes in the 
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Fig. 1. Drift coefficients b for the stochastic differential equation eq. (2.10) with n = 1 associated 
with the periodic solution (3.1a) for C > 1. 




Fig. 2. Drift coefficients b for the stochastic differential equation eq. (2.10) with n = 1 associated 
with the kink solution (3.1c), derived by setting C = 1. 

same neighborhood. The periodicities or localized structures affect the behavior of the sample 
paths of Xf. 

4. Numerical Results of the Stochastic Equations 

In this section, wc shall analyze eq. (2.10) numerically and study the time evolution of Xt, 
associated with the solution eq. (3.1). For all the simulations, we shall take m = 1. The order 
n, which determines the figure of the diffusion coefficient, is selected to be 1 and 5. We note 
that all the sample paths for each case start from a fixed point, although the initial points 
should be distributed in accordance with the initial distribution p(x,t = 0). However, as we 
will see shortly, the effect of this omission remains only in the early stage and disappears after 
the stochastic variable Xt is captured by the region where the stochastic effect is sufficiently 
strong. 
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4-1 Time evolution of the stochastic variable associated with the periodic solutions for n = 1 
In Fig. 5, five sample paths for tlie solution eq. (3.1a) are sliown togetfier. Tfie parameters 
in 6, a and (p are selected as 

C= 1.2, a = 1, ^ = 0.5, (4.1) 

and the initial condition is set as Xq = 0. The paths fluctuate around their initial points in 
the early stage. Although they split into some groups of paths, the variable Xt changes along 
several lines parallel to each other. The average velocity of the motion of Xt is observed to be 
approximately the same as that of the periodic solution. 

4-2 Time evolution of the stochastic variable associated with the kink solution for n = 1 

Another example is a stochastic process associated with the 1-kink solution eq. (3.1c). 
The coefficients are given by 

am , m(x — vt) , , , 

b{x,t)=v ;^tanh ^ (4.2) 
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Fig. 4. Examples of the coefficients of the stochastic differential equation eq. (2.10) with n = 5, 
associated with the kink solution (C = 1). (a) Drift coefficient b. (b) Diffusion coefficient a. 

We start with a set of parameters, 

C = 1, a = 1, v = -0.5, (4.3) 

and the initial condition Xq = —3. The results of the five paths are shown in Fig. 6. As we 
can see from the figure, all the paths show similar features. The paths start from x = Xq, 
linger around the initial point owing to the effect of Brownian motion, and finally travel along 
a line on the X-t plane. Although some of the paths display large deviations, or stay around 
simple values occasionally, the average incline of each path coincides with the velocity of the 
peaks of the periodic solution of 

4-3 Sample paths in the case associated with the kink solution for n = 5 

At the end of this section, we shall show the behavior of Xt for n = 5. Because the explicit 
forms of the drift coefficient for periodic solutions in Fig. 3(a) are similar to the iteration of b 
associated with the kink solution in Fig. 4(a), we chose the same parameters listed in eq. (4.3), 
with the initial condition Xq = —1.5. The behaviors of Xf are shown in Fig. 7. The paths 
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Fig. 6. Sample paths for Xt associated with the kink solution of the SG equation with n = 1. 



initially draw curves near a simple line with small fluctuations, and some of the paths undergo 
nearly the same amount of position shifts and begin oscillating with rather large amplitudes. 
In both of the initial and final states, the average incline of the paths arc observed to be the 
same. Although each path of Xt fluctuates intensively, and shifts from the initial line, its mean 
path is parallel to a set of lines in the initial state. The mean paths in the final state trace 
the line drawn by the center of the kink that passes through the origin, while the paths in the 
initial state have a common intercept corresponding to Xq. The time when a path falls into a 
state with a strong oscillation differs from time to time. 

5. Asymptotic Behaviors of Sample Paths 

In the numerical results shown in the previous section, the behaviors of the sample paths 
for Xt in sufficiently large time regions imply the existence of some deterministic process 
that drives the mean movements of the stochastic variable, and we shall roughly describe this 
effect in terms of the equation of motion generated from eq. (2.10). When we omit the second 
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Fig. 7. Sample paths for Xt associated with the kink solution of the SG equation with n = 5. 

term in eq. (2.1), we have a nonstochastic time evolution law of Xt, and we can find that the 
relation dXt = b{Xt,t) dt holds. This means that the variable Xt represents the position of a 
particle whose velocity is b{Xt,t). Since we consider traveling wave solutions of eq. (1.1) and 
the drift coefficient given by eq. (2.9), the dependences of b{x,t) on space and time variables 
have the same structure as those of (p, which means b = b{x — vt). Differentiating the relation 
X = b{x — vt) , we have the relation 

X = ^b{x, t) = b' ■{x-v) = b'x - vb'. (5.1) 

When we consider a small region around a given point of x, and approximate the graph of b 
by its tangent whose incline is /3, eq. (5.1) reduces to 

x = (3x- /3v. (5.2) 

This equation is the same as that for a particle with resistivity proportional to its velocity 
under a constant force. For negative /3, the motion results in a state of constant velocity v, 
whereas the value of x diverges for positive /?. Thus, there appear attractors for x at points 
where both b{x — vt) = v and b'{x — vt) < are satisfied. 

As an example of this explanation, let us first roughly describe the case associated with the 
kink solution forn = 1. The exact figure of the drift coefficient is given in Fig. 2. The coefficient 
b is almost constant in regions far from the center of the kink, whereas it conspicuously 
decreases with x around the center. In far regions from the center of the kink, the variable x 
moves almost at a constant velocity toward the attractive region. After x is caught in the area 
where the time evolution of the variable x is approximated by eq. (5.2), the path approaches 
the trajectory of the center of the kink, since the attractive region also travels at a velocity v. 

Next, we shall explain the behaviors of the paths for each case of the simulation in detail. 
In the case of n = 1 and C > 1, there appear periodic intervals, whose centers are the attr actor 
of X. The initial point, Xq = 0, is one of the repulsing points, and Xt travels away from the 
initial point at an average velocity of b{Xt,t). The motion of Xt shows a small fluctuation 
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due to the term adBf. Finally, each path Xt is trapped in one of the attractors and moves 
at an average velocity v. The three paths in Fig. 5 are located around the trajectories of the 
attractors. Since the initial velocity of Xt has a negative value, b{XQ, 0) = v < 0, the number 
of paths in the region of large Xt tends to be smaller than that for smaller Xt- 

The other case of n = 1, the result associated with the kink solution, can be explained 
in the same way. Only a single attraetor appears in this case, and all the paths fall into the 
line X = vt in Fig. 6. At the initial position Xq = —3, b is approximately b ~ 0.2, and Xt 
increases gradually with stochastic fluctuations. After the point t ~ 4, the paths are caught 
by the attraetor near the center of the kink, and Xt is driven by the attraetor. 

In the case n = 5 and C = 1, we should consider the property of the diffusion coefficient 
as well as the drift coefficient, because a is not a constant. As in the ease of n = 1, the point 
b{x — vt) = with 6' < plays the role as an attraetor, but the intensity of the stochastic 
effect is valid only around the center of the kink. When the Xq is sufficiently far from the 
center, Xt shows no stochastic effect and travels at a velocity v, and Xt varies with little 
stochastic oscillation. In the case shown in Fig. 7, b is approximately 6 ~ ti = —0.5 and 
a is rather small but cannot be negligible near Xq = —1.5. Thus, Xt initially moves at a 
velocity v approximately, with little fluctuation. After Xt is caught by the attraetor, it shifts 
its trajectory near the center of the kink and begins to oscillate at a large amplitude because 
of large values of a. The point when Xt begins to oscillate is determined purely in a stochastic 
way, and this causes some of the paths to fluctuate extremely (Fig. 7), but others travel in a 
rather deterministic way. 

6. Discussions and Conclusions 

In this paper, we have presented new stochastic equations, associated with the traveling 
wave solution of the SG equation. By virtue of the structure of the unknown function of 
the SG equation, as a traveling wave, and the assumption that the density function p{x,t) is 
proportional to <^x, we have found that the coefficients of the stochastic equation eq. (2.1) can 
be expressed in terms of the derivatives of 0, because the coefficients in the equation for Xt 
should not have singularities. We have also performed numerical analyses of the time evolution 
of Xt under some initial conditions. The paths drawn by Xt show characteristic appearances, 
namely, the average Xt asymptotically coincides with the velocity of the traveling wave. These 
behaviors of the paths can be explained by considering attractive points generated from the 
drift coefficient, which moves at the same velocity as the solution itself. Thus, the paths Xt 
are observed to be driven by the solution of the SG equation. 

The results derived in this paper are distinctive in that we can consider stochastic motions 
of a particle, where the density function is expressed in terms of the traveling wave solution 
of a solvable equation. Since the traveling wave solutions of the SG equation are well known, 
this feature enables us to derive various physical quantities related to stochastic evolutions 
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accurately. In addition, we can expect several applications of the result presented in this 
paper. The behaviors of solitons under external effects such as impurities"*^^' or fluctuations 
in metamaterials^^ have attracted much interest and have been studied in various systems. 
Although many results have been derived from these studies, these studies focus on effects that 
act on solitons; the stochastic effect caused by solitons has not been studied so far. Applying 
the results in this paper to the analysis of phenomena in real materials, we can find how 
solitons affect original systems as feedbacks, and how particles excited in media travel while 
interacting with solitons. 

Since the SG equation is related to many physical systems, the derived model can provide 
good descriptions of such systems including stochastic effects. For example, the mechanical 
model of the motion of connected rotators, which are applied forces due to twisted connection 
proportional to twist angle under gravity,^ or a dislocation in sinusoidal periodic potentials"'^'' 
with heat fluctuation, and the optical waves in two-level atoms^"*^ with heat excitations, are 
such candidates. Therefore, we can expect that the stochastic processes derived in this paper 
can be observed experimentally. Such experiments in actual media or derivations of stochastic 
equations from physical models are challenging problems. 

At the end of the paper, we would like to discuss extensions of the derived results. We 
used the special structure of traveling waves to derive the stochastic differential equation, 
namely, their coefficients b and a in explicit forms. Because of this restriction, the result is 
not applicable to the case associated with a solution that has a time-dependent profile such 
as a two-soliton solution. When we construct stochastic processes for general solutions, we 
can consider stochastic phenomena associated with the SG equation under various conditions. 
Another extension is to consider stochastic equations for other integrable equations. For exam- 
ple, the nonlinear Schrodinger (NLS) equation is one of the well-known models in integrable 
systems. The method of stochastic quantization is known to be effective way for connecting 
the Schrodinger equation with a stochastic differential equation, and the stochastic equation 
associated with the NLS equation can be used to describe the motion of particles in various 
physical models such as plasma physics or nonlinear optics. There are many applications to 
be considered, and these are future problems to be examined. 
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